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Abstract 

The effects of magnetism on the Bain transformation of a-phase FeNi systems are investigated 
by using the full potential linearized augmented plane wave (FLAPW) method based on the gen- 
eralized gradient approximation (GGA). We found that Ni impurity in bcc Fe increases the lattice 
constant in ferromagnetic (FM) states, but not in the nonmagnetic (NM) states. The shear mod- 
ulus G and Young's modulus E of bcc Fe are also increased by raising the concentration of nickel. 
All the compositions considered show high shear anisotropy and the ratio of the bulk to shear mod- 
ulus is greater than 1.75 implying ductility. The mean sound velocities in the [100] directions are 
greater than in the [110] directions. The Bain transformation, which is a component of martensitic 
transformation, has also been studied to reveal that Ni^Fei-x alloys are elastically unstable in the 
NM states, but not so in the FM states. The electronic structures explain these results in terms 
of the density of states at the Fermi level. It is evident that magnetism cannot be neglected when 
dealing with the Bain transformation in iron and its alloys. 
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I. INTRODUCTION 



Nickel and iron form substantial solid solutions over the complete range of compositions. 
The addition of nickel to iron enhances the strength and toughness of iron. In its a-phase, 
the iron-nickel alloy forms the base of structural steel; henceforth, the mechanical properties 
of these alloys received more attention as compared to other properties. The body-centered 
cubic (bcc) a-phase is stable from pure iron to ~ 10 % Ni and increasing the Ni concentration 
in a-Fe reduces its stability with respect to 7-Fe which is why nickel is designated a 7- 
stabilizer.- Fe-Ni alloys are also of interest in connection with Invar effect and martensitic 
transformation at low temperatures.-"^ Fe-Ni alloys are amongst the most studied magnetic 
materials and are important in understanding the mechanical and magnetic properties of 
steels. Considerable efforts have been made to understand the stability of magnetic materials 
in terms of magnetism and Bain transformation," but with the focus on Bain path of either 
pure bcc Fe, or fee Fe, or hep Fe, rather than on the consequences of the presence of nickel 
as a solute. 

Much attention has also been given in the past to Fe-Ni alloys near Invar composi- 
tions or ordered compounds (FesNi, FeNi, or FeNis), but not dilute alloys with less than 
10 % nickel .-i^>^ii^ Recently,— we studied the electronic and magnetic structures of transition 
metal impurities in bcc Fe where we also found that Ni impurity enhances the magnetism 
of bcc Fe. To further elucidate the effect of Ni on iron, we examine here the elastic proper- 
ties, thermodynamics, and the effects of magnetism on the Bain path using first-principles 
calculations. 

II. COMPUTATION MODEL AND METHODS 
A. Computation Model 

We considered a 3 x 3 x 3 supercell of the primitive cell of bcc Fe, which contains 27 
atomic sites. The centered Fe atom is replaced by the Ni atom; we denote the model as 
NiiFe26- If we replace the central Fe atom by Ni, then the space group is ImSm (space group 
#229), and the Ni atom is at the (2a) Wyckoff position, its 8 nearest neighbor Fe atoms 
(denoted as Fel) are on the (16/) Wyckoff sites, the 6 second nearest neighbors (denoted as 
Fe2) are on the (12e) Wyckoff sites, and the 12 third nearest neighbors (denoted as Fe3) are 
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on the (24/i) Wyckoff sites. A 2 x 2 x 2 supercell of a conventional cubic cell of bcc Fe, which 
contains 16 atoms, was also considered and this system is denoted as NiiFeis. The supercell 
is a simple cubic unit cell, with lattice constant twice that of the bcc structure. The space 
group is Pm3m (space group ^^221), and the Ni atom is on the (la) Wyckoff position. The 
neighboring atoms are on {8g), (3c), {3d), and (16) sites, in order of increasing distance from 
the central Ni atom. The corresponding Ni concentrations are 3.7 and 6.25 at.% for 3 x 3 x 3 
and 2x2x2 supercells, respectively. Pure bcc Fe, bcc Ni, and fee Ni as references were 
also considered. 



B. Electronic structure calculation method 



The Kohn-Sham equation was solved self-consistently in terms of the total energy all- 
electron full-potential linearized augmented plane-wave (FLAPW) methodi^*^ based on the 
generalized gradient approximation (GGA) for the exchange-correlation potential.— The 
integrations over the three dimensional Brillouin zone (3D-BZ) were performed by the im- 
proved tetrahedron method^ over a 13 x 13 x 13 Monkhorst-Pack mesh^ in the 3D-BZ, which 
corresponds to 84 k-points inside the irreducible wedge of the 3D-BZ. The linearized aug- 
mented plane-wave (LAPW) basis set was expanded using a plane wave with an energy 
cutoff at 4 (27r/a), where a is the lattice parameter. Lattice harmonics with / < 8 were 
employed to expand the charge density, potential, and wave functions inside each muffin-tin 
(MT) sphere with the radius of 2.20 a.u. for both Fe and Ni atoms. The star-function 
cutoff of 16.73 {2ti /a) was employed for depicting the charge density and potential in the 
interstitial region. The core electrons were treated fully relativistically, and the valence 
electrons were treated scalar relativistically. To ensure the orthogonality between core and 
valence states, we employed the explicit orthogonalization (XO) method.— Self-consistency 
was assumed when the difference between input and output charge (spin) density is less than 
1.0 X 10~^ electrons/a. u.^ The convergence of these computational parameters was carefully 
checked.— 
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C. Methods for elastic properties 



Calculations were carried out on nonmagnetic (NM) and ferromagnetic (FM) states, to 
find the optimized lattice constant by fitting the total energy data at various lattice constants 
to the Birch-Murnaghan equation of state (EOS),— which also provides the bulk modulus 
B. 

Since the considered systems are cubic, there are only three independent elastic constants; 
Cii, Ci2, and C44. These elastic constants can be determined by calculating the total energy 
as a function of the shears described below. To determine Cn and C12, we considered a 
volume-conserved tetragonal strain in such a way as to modify the cubic crystal axes by 
applying the following strain matrix, 

^ 1 + 5 ^ 

1 + 5 
^ (1 + 5)^2^ 

where 5 is the amount of strain ±0.03 imposed on the crystal. The change in the strain 
energy density {u) as a function of strain is given by 

n = 6CV + 0(<5^), (1) 

where C is the tetragonal shear constant defined by (Cii-Ci2)/2. By calculating C and 
bulk modulus relation B =|(Cii+2Ci2) one can estimate Cu and C12. We must keep in 
mind that the bulk modulus B was calculated by using the EOS.— To determine C44, the 
following volume-conserved orthorhombic distortion 

lis \ 

5 1 
^0 (1-5^)-^^ 

was considered and C44 can be calculated from 

u = 2Cii5^ + 0{5^). (2) 

The above elastic constants can be put into a more general way, i. e, the single-crystal shear 
moduli for the {100} plane along the [010] direction and for the {110} plane along the [llO] 
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direction are simply given by Gjioo} = ^"44 and Gjno} = {Cn — Ci2)/2, respectively.— Once 
C12 and C44 are known, it becomes possible to calculate the Cauchy pressure Cp — C12 — C44. 

The shear modulus G, Young's modulus E, Poisson's ratio u, and shear anisotropy factor 
A for polycrystalline aggregates can be calculated from the elastic constants. The shear 
modulus is bounded by the Reuss' Gr modulus and Voigt's Gy one which represent their 
lower and upper limits, respectively.— For cubic lattices these moduli are determined by 

^ _ 5(Cii — Ci2)C44 . . 

''"4C44 + 3(Cn-Ci2)' ^ ^ 



and 



Cy = ~ + ^^^^^ . (4) 



Despite wide usage, neither Russ' nor Voigt's relation is believed to be exact. Hill2^ suggested 
an averaging by arithmetic mean of Gr and Gy 

G=^iGn + Gv). (5) 
The bulk and shear moduli are used to calculate the Young's modulus 

E = ^ (6) 

3B + G' ^ ^ 



and the Poisson's ratio 

3B-E 



6B 

The shear anisotropic factor A can be calculated by using the following relation^i 



(7) 



A = ^=r^. (8) 



D. Methods for thermodynamic data 

Using the equilibrium lattice constant, the formation enthalpy AH per atom of Ni^Fe^ 
was calculated as follows: 

AH = ■ , (9j 

m + n 

where if(Ni„Fem) is the enthalpy of Ni^Fcm with m = 26, 15 and n = 1, and H{Fe) and 
H(Ni) are the total energies per atom of the ground states of bcc Fe and fee Ni, respectively. 



The Debye modeP^ was used to calculate the mean isotropic velocity which is given 



1 -1/3 



by 



where the longitudinal w_l and transverse v\\ velocities are given by 

= V(5 + 4G/3)/p, 

and 



(10) 



(11) 



(12) 



where p denotes the density of material. 

Once we know all the elastic constants of a cubic system, we can also examine the 
behavior of sound velocities in different crystallographic directions, e.g., [100], [110], or [111] 
directions.— The velocity of longitudinal [vi) and transverse {vt) elastic wave in the [100] 
direction is given by 



Vl = I — 1 , Vt 
Similarly, for the [110] direction 

Cii + Ci2 + 2C44^ ^ 

Vl 



P 



Vt 



2p y ' V 2p 

and for the [111] direction, the velocities can be expressed as 



(13) 



(14) 



Vl 



Cu + 2Ci2 + 4a 



-■44 



1/2 



Vt 



Cn — C12 + C44 



1/2 



3p J \ 3p 

Vl and Vt were used to distinguish it from the isotropic velocities given in Eqs. f|TQ] - [T2| 



(15) 



III. RESULTS AND DISCUSSION 



A. Formation Enthalpy 

The calculated lattice parameters of bcc Fe, Ni, fee Ni, bcc NiiFe26, and NiiFcis in the 
NM and FM are presented in Table [11 The lattice constant of bcc Fe is calculated to be 
2.76 A (2.83 A) in the NM (FM) state which is comparable with the previous calculations of 
2.84 A in the FM state and 2.76 A in the NM state^^ and the experimental observations of 
2.87 A.— The equilibrium lattice constant of NiiFe26 is determined to be 2.76 (2.84) A in the 
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TABLE I. The calculated lattice constant (a) in units of A of bcc Fe, bcc Ni, fee Ni, and bcc 
NiiFe26 and NiiFeis. Both the nonmagnetic (NM) and ferromagnetic (FM) results are listed. The 
experimental values (Expt.) of the bcc Fe,— bcc Ni,— fee Ni,— and bcc NiiFe26,— and the other 
theoretical results^i (Other) are also listed for comparison. 





bcc Fe 


bcc Ni 


fee Ni 


NiiFe26 


NiiFeis 


NM 


2.76 


2.79 


3.51 


2.76 


2.76 


FM 


2.83 


2.80 


3.52 


2.84 


2.85 


Expt. 


2.87 


2.82 


3.52 






Other 


2.83 


2.80 


3.52 







NM (FM) state. The corresponding values for NiiFeis are determined to be 2.76 (2.85) A 
in the NM (FM) state. Table [T] shows that the lattice constant of bcc Fe increases with the 
Ni concentrations. This is consistent with experimental observations, which show that Ni 
as a solute slightly expands the lattice constant of bcc Fe.— It is interesting to find that the 
lattice parameter of bcc Fe is not affected by the Ni addition in the NM states; this can be 
attributed to the magnetovolume effect.— Therefore, Ni expands the lattice of bcc Fe in the 
FM state, but not in the NM state which suggests that the ferromagnetic interactions are 
responsible for the increment of the lattice constant of Fe with the Ni concentrations. The 
present calculated values are slightly smaller than those determined at room temperature 
as might be expected from thermal expansion,— but it is encouraging that the trend of the 
lattice parameter as a function of the nickel concentration is correctly reproduced.— It is 
also noticeable that ab initio calculations usually give trends and it may underestimate or 
overestimate the lattice parameters. Our previous calculations on bcc Fe-based materials 
show that atomic relaxation around the impurities is negligible.—*^ 

Using the optimized lattice parameters, AH was calculated using Eq. iQ and the results 
are given in Fig. [H which shows AH of fee NiFe systems taken from Ref. 
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Note that 



AH of a system is nothing more than the total energy of the system at zero pressure and 
zero kelvin at the corresponding equilibrium lattice parameter. The formation enthalpy 
is measured relative to bcc Fe and fee Ni. We see that the small addition of Ni in bcc 
Fe increases AH. Based on previous^^ and present study,it can be summarized that the 
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FIG. 1. Calculated formation energy (kJ/mol atom) of bcc Fe, Ni, NiiFeis, NiiFe26, and fee Ni 
versus the Ni concentration. Open symbols represent our work (this) and filled symbols show other 
calculated values (other) taken from Ref. liol . Circles (triangles) indicate the bcc (fee) systems. 



addition of Ni to Fe stabilizes fee struetures and destabilize the bcc structures, consistent 
with the experimental phase diagram.- 



B. Elastic Properties and Thermodynamics 

In order to obtain the elastic properties, ±0.03 strains were imposed as discussed in 
Eqs. (II]) and (E]). The calculated bulk moduli and elastic constants of bcc Ni^^Fei-a; and fee 
Ni are shown in Table [Tll The calculated elastic constants of bcc Fe and fee Ni are close to 
the experimental values, and confirm that the number of k-points and the number of basis 
functions, used in these calculations, are sufficient to reproduce the experimental data.— All 
the systems have positive elastic constants which satisfy the mechanical stability condition 
of a crystal.— However, bcc Ni is mechanically unstable due to the negative estimated value 
of Gjioo}- Although it is mechanically unstable, bcc Ni has been successfully achieved as a 
thin film on GaAs substrate.— If Gjno} is larger than Gjioo}? then it is easier to shear on 
the {110} plane than on the {100} plane. It is noticeable in Table HT] both the bulk modulus 
and elastic constants(Cii and C12) of bcc NiFe26 are increased. With the addition of Ni in 
bcc Fe, Giiio} is increasing while Gjioo} is decreasing. 

Table UTTIlists the shear G and Young's modulus E, B/G, Cauchy pressure Gp, anisotropic 

8 



TABLE II. The calculated elastic constants and bulk moduli in the units of GPa of bcc Fe, fee Ni, 
bcc Ni, and bcc NiiFe26- The experimental (Expt.) values for the bcc Fe^^ and fee Ni^^ are given 
for comparison. 





Cii 


Cl2 


C44 


Cjioo} 


B 


bcc Fe 


257 


134 


95 


61.00 


175.00 


Expt. 


242 


136 


112 


53.00 


171.33 


fee Ni 


276 


156 


130 


60.00 


196.00 


Expt. 


248 


155 


124 


46.50 


186.00 


bcc Ni 


141 


211 


154 


-35.00 


187.67 


NiiFe26 


267 


149 


102 


59.00 


188.33 


NiiFeis 


243 


127 


110 


58.00 


165.67 



TABLE III. The calculated shear modulus G (in units of GPa), Young's modulus E (in units of 
GPa), B/G^ Cauchy pressure Cp (in units of GPa), anisotropy factor A, and Poisson's ratio v of 
bcc Fe, bcc Ni, fee Ni, and bcc NiiFe26 and NiiFeis. 





G 


E 


B/G 


Cp 


A 


V 


bcc Fe 


79.80 


207.82 


2.19 


39 


1.54 


0.30 


bcc Ni 


-27.18 


-85.67 


-6.90 


57 


-4.40 


0.58 


fee Ni 


95.32 


246.07 


2.06 


26 


2.17 


0.29 


NiiFe26 


81.89 


214.57 


2.30 


47 


1.73 


0.30 


NiiFeis 


85.08 


217.94 


1.95 


17 


1.90 


0.28 



factors A, and Poisson's ratio v for Ni^Fei.^.. Note that G is bounded either by Gy or Gji 
and, therefore, G is in an approximated average sense. It is a common practice to consider 
G as an indication of mechanical properties of materials.—"^ It is known that the hardness 
and strength of materials are related to their elastic moduli, such as -E, -B, and G.— The 
bulk modulus of bcc Fe is slightly increased by 3.7 at. % Ni. It is noticeable that both E 
and G increase with Ni. On the other hand the negative values of bcc Ni demonstrate that 

9 



TABLE IV. The calculated sound velocities (in units of km/s) and densities (in units of kg/m'^) of 
bcc Fe, fee Ni, and bcc NiiFe26 and NiiFeis. 





v± 


^11 




P 


bcc Fe 


5.86 


3.12 


3.48 


8190 


fee Ni 


6.01 


3.27 


3.64 


8940 


NiiFe26 


6.05 


3.18 


3.55 


8120 


NiiFeis 


5.89 


3.25 


3.62 


8040 



this phase violates the mechanical stability condition.— All these materials have high shear 
anisotropy factor and anisotropic behavior as increasing Ni concentrations. Pettifor— has 
suggested that the Cauchy pressure could be used to describe the angular character of atomic 
bonding in metals and compounds. If the bonding is more metallic, the Cauchy pressure 
will be positive. It is believed that ductile materials (such as Ni or Al) have positive values 
and brittle materials (such as Si) have negative values of Cauchy pressure.— This can be 
seen in Table IIIII where Nij.Fei_a; have positive values and are ductile, consistent with the 
metallic behavior discussed below. 

To further shed light on whether Ni^^Fei-a; are brittle or ductile, we used simple Pugh 
relations that link empirically the plastic properties of metals with elastic moduli by B/GM 
If the B/G ratio is greater than 1.75, the material behaves in a ductile manner, otherwise 
in a brittle manner, as demonstrated by first-principles calculations.—"^ It is seen that the 
B/G ratio is greater than 1.75 for bcc Nia;Fei_a; pointing that all our considered systems are 
ductile as expected. The other factor that measures the stability of a crystal against shear 
is Poisson's ratio u, and it has been shown^!^ that brittle materials have higher Poisson's 
ration such as NiAl (z/ = 0.41), while the ductile materials have Poisson's ratio u ~ 0.30. 
Now it is clear from Table UTTI that the Poisson's ratios of bcc Ni^-Fei-^ are ~ 0.30. Therefore, 
bcc Nia;Fei_2; enter into the class of ductile materials. At low temperature the vibrational 
excitations arise solely from acoustic vibrations and the calculate isotropic mean velocities, 
using Eq. (ITOl) , are given in Table HVl 

We emphasise that the modulus ratio is just one of many indicators of the propensity 
of a material to be brittle;— crystal defects, the weakening of grain boundaries by chemical 

10 



TABLE V. The calculated sound velocities (in units of km/s) of bcc Fe, fee Ni, and bee NiiFe26 
and NiiFeis in different directions. The values of Vm are calculated by using Eq. (fTOj) . 



[100] [110] [111] [100] [110] [111] [100] [110] [111] 



bcc Fe 5.60 5.96 6.07 

fee Ni 5.56 6.22 6.43 

NiiFe26 5.73 6.18 6.32 

NiiFeis 5.50 6.06 6.23 



3.41 2.74 2.98 

3.81 2.59 3.05 

3.54 2.70 3.01 

3.70 2.69 3.06 



4.34 3.72 4.00 

4.66 3.57 4.13 

4.49 3.69 4.06 

4.55 3.67 4.11 



segregation and many other factors can play a seminal role in determining the macroscopic 
properties. We found that the calculated elastic constants are in agreement with the available 
data, and the knowledge of elastic constants permits the sound velocities to be estimated as 
a function of the direction. Simply using Eq. ( IT314T51) . the sound velocities of all the systems 
considered here were calculated, and the velocities of bcc Nia;Fei_a; in the [100], [110], and 
[111] directions are given in Table IVl 

One can see that Vm in the [100] direction, which is also an easy axis of magnetization 
of bcc Fe, is increasing with the Ni concentrations. On the other hand, Ni concentrations 
decrease fm in the [110] direction. It is clear to see that has larger values in the [100] di- 
rection than in the [110] direction. The different values of in the different crystallographic 
directions confirm the anisotropic behavior of bcc Nia.Fei_a;. 

Before going to discuss the Bain path, it is emphasized that the elastic constants calcula- 
tions did not allow internal coordinate or stress-tensor optimization at each lattice distortion. 
Since lattice distortions were small ~ ±3%, this constraint may preserve the stable states 
and the effect of atomic relaxation (in either direction, i.e., x,y,z) is found to very small. 
However, we also performed some test calculations on the internal coordinate relaxations 
in the distorted systems, but they were hardly distinguishable with the unrelaxed ones. It 
is believed that atomic relaxation will not have much effect (qualitatively) on the elastic 
properties and such relaxation will not change the main conclusions drawn in this paper. 
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c/a c/a c/a 

FIG. 2. Total energy (in eV units) as a function of c/a for (a) bcc Fe (b) NiiFe26, and (c) NiiFeis. 
Filled (Open) circles represent the total energy in NM (FM) states. 



C. Magnetism on the Bain deformation 

Let us examine the effects of magnetism on the martensitic transformation in terms of 
the Bain deformation.— One should not be confused with the reversibihty of martensitic 
transformation by crystal symmetry,— because we ignored the effects of entropy totally, 
which causes the irreversible processes.— The mechanical stability of the cubic phase under 
tetragonal distortion was evaluated by calculating the total energy as a function of volume 
conserving tetragonal strain at the optimized lattice parameters of bcc Nia;Fei_a;. The total 
energies in the NM and FM states as a function of c/a ratio at constant volume are presented 
in Fig. [2j It is seen that in the NM states the total energy has a negative curvature around 
c/a = 1.0, which is a signature of mechanical instability.— It is further expected that the 
structure instability will be accompanied by a softening of phonons.-i^i^i^ These figures 
clearly show that when the cubic symmetry of the bcc phase is broken by the tetragonal 
distortion in the NM state, the total energy is decreasing and indicating a phase transition to 
the other phase; fee in this case. The total energy is minimized around c/a = a/2 (fee phase) 
and become mechanically stable. This trend can be seen in all these bcc FeNi systems. 

However, the situation was drastically changed when we repeated the same calculations 
in the FM states. In the FM states, these systems are mechanically stable in the bcc phases 
and unstable in the fee phases. It is possible on the basis of these observations to conclude 
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that the bcc phases are stable in the FM states, and fee phases in the NM states. The bcc- 
fcc phase transition is prohibited by the magnetic contribution to the total energy. If these 
systems were not magnetic, then one could expect bcc-fcc phase transitions in these bcc 
FeNi systems. This further points out the role of magnetism in the phase stabilities of FeNi 
systems, i.e., the bcc phases are stabilized by ferromagnetism. These findings show that 
magnetism contributes to both structural and elastic stabilities of materials. Nevertheless, 
magnetism can either stabilize an elastically unstable material, as discussed above, or cause 
a stable structure to become elastically unstable, e.g., hep Fe. A similar Bain transformation 
path was also observed for bcc-fcc Fe where it has been shown that the local minima on the 
Bain path was changed by changing the lattice volume.-*^ Therefore, the Bain path can be 
altered slightly by changing the lattice volume. 

It is worth to mention how our magnetic model describes the reality of alloy system. 
Our comparison of the ferromagnetic and nonmagnetic solutions is simply to emphasise the 
stability of the former state. The comparison should ideally be between the ferromagnetic 
and paramagnetic states which we cannot easily access. Alloys differ from the corresponding 
chemical compounds by introducing the concept of atomic thermal aggitation,— which can be 
interpreted as an Ising ferromagnet.— Of a given atomic configuration at a finite temperature, 
the Ni impurity modifies the the magnetic moments of neighboring Fe atoms. In addition, 
Fe-Ni system is a representative itinerant ferromagnet^ described by the Stoner model.— 
A natural long-range oscillations, so-called the Ruderman-Kittel-Kasuya-Yosida (RKKY)^ 
and Friedel^ oscillations, may occur. However, the thermal aggitation will change the 
atomic configurations dynamically and it will diminish the effects of the RKKY and Friedel 
oscillations, by superposition of the random phases of those oscillations. Hence, we can 
expect that the averaged magnetic moments of atoms will not much differ from those of a 
representative atomic configuration. Although there are good methods for simulating the 
magnetic moments of disordered states,— the problem of alloy magnetism is still open. In 
addition, bcc iron matrix screens well such oscillations of the magnetic moments within the 
nearest neighbor distance.— "^"^ Therefore, our model should not change the essence of the 
conclusions, since the ferromagnetic state of iron and its dilute alloys is well-established to 
be so much more stable than the paramagnetic state. 

13 




1.5 I ' ' ' ' ' ' ' ' 1 

0.92 1 1.08 1.16 1.24 1.32 14 1.48 1.56 

c/a 

FIG. 3. Total density of states (DOS) per atom at Ep of bcc NiiFeis and NiiFe26 as a function of 
c/a. Filled (Open) squares represent the DOS of NiiFeis (NiiFe26) in the NM states. 

D. Electronic structures on the Bain deformation 

The observed stabihty/instabihty in fcc/bcc phase ultimately have an underlying elec- 
tronic origin. The electronic structures of these bcc systems and the calculated total elec- 
tronic density of states (DOS) per atom at the Fermi level [Ep), denoted as nlEp), in the 
NM state, are given in Fig. [Hi The bcc phase, Nij.Fei_a; has a minimum n{Ep), and as 
it is deformed to body-centered tetragonal (bet) then n (Ep) suddenly increases and has a 
maximum value around c/a = 1.08. The lowest n{Ep) can be seen near c/a = This 
shows that the bcc structure is mechanically unstable due to large n (Ep). 

To reveal the effect of tetragonal distortion on the local DOS, the atomic projected local 
DOS of NiiFe26 and NiiFeis is shown in Fig. H] and Fig. |5l respectively, which show the local 
DOS for c/a = 1.0 and c/a = ^/2. The DOS shows metallic behavior which also follows 
the Pettifor— suggestions about the metallic character and positive values of the Cauchy 
pressure. One can clearly see the reduction (~ 50%) of n (Ep) for the c/a = \/2 cases. In 
both NiiFe26 and NiiFeis, n{Ep) is dominated by the Cg orbitals for c/a = 1.0, however, 
the orbitals contribute to n (Ep) for c/a = It is also clear from these figures that 
Ep of cubic Nij-Fei-a; is located on a peak of high DOS and the tetragonal distortion is 
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FIG. 4. (Color online) Calculated spin-unpolarized impurity-atom-projected local density of states 
of NiiFe26 for (a) c/a = 1.0 and (b) c/a = \/2. Solid (dotted) lines represent the t2g (eg) states, 
whereas the bold solid lines show the total DOS per atom. The Fermi energy (Ep) is set to zero. 



leading to a decrease of n [Ep) in the tetragonal phase. The strong reduction in n [Ep) is an 
indicative of a states shifts from the Fermi surface, i.e., an instability associated with the 
peak in n (Ep). Note that in the absence and presence of distortion, the value of n {Ep) for 
bcc NijjFei-a; is within the Stoner limit.— In metals, in general, a lower n (Ep) corresponds 
to lower electronic kinetic energy. This reduction gives a substantial negative contribution 
to the electronic energy and consequently the total energy decreases upon lattice distortion 
as shown above, yielding a negative curvature of the shear modulus for bcc FeNi systems. 
This mechanism is responsible for the elastic instability of bcc Nij-Fei-a, in the NM states. 
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FIG. 5. (Color online) Calculated spin-unpolarized impurity-atom-projected local density of states 
of NiiFeis for (a) c/a = 1.0 and (b) c/a = \/2. Solid (dotted) lines represent the t2g {eg) states, 
whereas the bold solid lines show the total DOS per atom. The Fermi energy {E-p) is set to zero. 



IV. SUMMARY 

The elastic and thermodynamics of Ni impurities in bcc Fe is studied using first-principles 
calculations which showed that Ni impurities do not expand the lattice constant of bcc Fe 
in the nonmagnetic state, in contrast to the the ferromagnetic state. Nickel impurities 
improved the elastic properties of bcc Fe. The sound velocities of elastic waves were shown 
to be increasing in the [100] directions, and decreasing in the [110] directions with the Ni 
concentrations. We also investigated the Bain path of bcc Nia,Fei_2;, and it was observed that 
Ni^iFei-a; systems are elastically unstable for c/a = 1.0 in the nonmagnetic sates, whereas 
Nij-Fei-j; systems are stable for for c/a = a/2. The elastic stability was explained in terms 
of electronic structures of Ni^-Fei-j;. 
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